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Abstract 

Different mechanisms which may possibly explain existence of magnetic fields 
on astronomically large scales are described. A recently suggested model of 
magnetic fields generation slightly before hydrogen recombination is discussed 
in more detail. 



1 Introduction 



It is known from observations, mostly from Faraday effect, that there are mag- 
netic fields in galaxies with magnitude of a few micro-Gauss and coherence 
length of the order of galactic size, l ga i ~ ( a few) kpc, The energy density of 
such fields is close to the energy density of the cosmic microwave background 
radiation (CMBR): 

Pb = B 2 /8tt ~ p. ( ~ 10~Vm « 3 ■ 1(T 34 g/cm 3 (1) 

where p 1 is the energy density of CMBR and p m w 1.5 keV/cm 3 is the total 
cosmological mass/energy density. Though the magnitude of such fields is 
small in comparison with the fields of stars or even planets (see below) it is 
very difficult to understand their huge coherence length. Even more puzzling 
is a possible existence of intergalactic magnetic fields which are 3 orders of 
magnitude weaker but coherent at a Megaparsec scale, for a review see ref. © . 

For comparison, magnetic field of the earth is = 0.5 G, magnetic fields 
of solar type stars can be as large as 10 3 G, magnetic fields of white dwarfs may 
reach 10 9 G, and absolute champions are neutron stars with magnetic fields at 
the level of 10 13 G. 

Possible creation mechanisms of large scale galactic or intergalactic fields 

can be roughly separated into four classes (for reviews see ref. 

1. Conventional (astrophysical) mechanism based on stellar ejecta of mag- 
netic fields with subsequent magnetic lines reconnection (from different 
stars) to create a homogenous on large scale field component. This mech- 
anism was reviewed recently in ref. ^ 

2. Processes in the early universe which invoke inflation to stretch the char- 
acteristic scale of the field up to galactic or even larger scales. A compre- 
hensive list of references on different versions of these mechanisms can be 
found in the recent paper » . 

3. Phase transitions in the early universe during which strong magnetic fields 
could be created but at a very small scale. 

4. Relaxation of previously created inhomogeneities. In this process tur- 
bulent or laminar flow of primeval plasma with non-zero vorticity could 
be generated and due to different mobilities of charge carriers vortical 
electric currents producing magnetic fields would be created. Such mech- 
anism may operate either in the early or relatively late universe. 

All these mechanisms either generate magnetic fields at very small scales or of 
insufficiently large amplitude. In the first case "Brownian" type reconnection of 



magnetic field lines is necessary for creation of coherent magnetic fields at large 
distances. The amplitude of the held in the course of reconnection decreases 
as (hn/lfin) 3 where k n and i/j„ are the initial and final coherence lengths. 
Inflationary stretched fields are usually rather weak, though may have a very 

large coherence length. In both cases galactic dynamo® should amplify the 
seed magnetic fields up to the necessary magnitude. 



2 A comment on astrophysical mechanism 

The total energy of galactic magnetic field in a large galaxy, e.g. in the Milky 
Way, is equal to: 

£gal = ^RgalPB~WM Q (2) 

where Mq is the solar mass. Magnetic energy of a neutron star with radius 
R ns = 10 6 cm is equal to « Kr n M (B/lO 13 G) 2 . Magnetic energy of 
white dwarfs with R wd = 10 9 cm is £^ d « 1CT 10 M Q (B / '10 9 G) 2 . Thus about 
10 11 white dwarfs with magnetic field B — 10 10 G or 10 12 neutron stars with 
B = 10 13 G should be in the Galaxy to feed galactic magnetic field if the 
energy of the field were not lost in the process of line reconnection. Since 
the total number of stars in the Galaxy is about 2 ■ 10 11 it is hardly possible 
to meet this requirement. Of course this estimate is quite rough and more 
accurate considerations may be not so negative but still it shows that efficiency 
of stellar creation of galactic magnetic field should be very high to satisfy energy 
constraints. 



3 Generation of seed magnetic fields in the early universe 

It is well known that during inflation very long gravitational waves and large 
scale scalar field perturbations are generated. A natural question is: "why not 
electromagnetic fields?" The answer is that scalars and tensor fields are not 
conformally invariant even in the zero mass case, while photons are. Conformal 
invariance means that rescaling metric with an arbitrary factor b(t, r) and the 
fields by the same factor to a power depending on the spin of the field we arrive 
to formally the same action written in terms of new variables. According to the 

Parker theorem® conformally invariant fields are not generated in conformally 
flat space-time. Indeed it is known that cosmological Friedman-Robertson- 
Walker (FRW) background is conformally flat, i.e. after the proper coordinate 
choice the metric can be written in the form 



ds 2 =a(r,x) 2 (dr 2 -dx 2 ) , 



(3) 



Thus after conformal transformation with the factor b — a(r, x) to a proper 
power we can exclude FRW-gravity for conformally invariant fields. Fortu- 
nately, as it has been already mentioned, this cannot be done for scalar (in 
particular, inflaton) fields. Otherwise cosmological density perturbations would 
not be generated and we would not be here. 

Several possible ways to break conformal invariance of electrodynamics 
and to create seed magnetic fields have been discussed in the literature: 

1. New non- minimal interaction of electromagnetic field with gravity, pos- 
sibly not gauge invariant^ : 

C = CxRA^A" + C 2 R^A»A V + C 3 R^ af3 F^ F a(i + ... (4) 

where A^ is the electromagnetic vector-potential, F^ v is the Maxwell 
tensor, R^ a /3 is the Riemann tensor, R^ v is the Ricci tensor, and R is 
the curvature scalar. 

2. Interaction with a new hypothetical field, dilaton, 

C = -(l/4)e e F^F^ (5) 



3. Quantum conformal anomaly due to famous triangle diagram which leads 
to non-zero trace of the energy-momentum tensor of electromagnetic field 

and breaks in this way conformal invariance of electrodynamics » ; 

T£ = kF^F^ (6) 

where the constant coefficient re depends upon the rank of the gauge 
group and the number of charged fermions. For SU(N) with Nt number 
of charged fermions it is: k — (a/ir)(UN/3 — 2Nf/3). 

In all these models electromagnetic waves could be generated at inflation- 
ary stage and sufficiently large magnetic fields at very large scales would be 
created to serve as seeds for galactic fields, if one takes appropriate values of 
of coupling constants in the Lagrangians or sufficiently large number of charge 
particles in the third case. Unsatisfactory features of this approach are intro- 
duction of new fields or interactions, though the second case looks quite natural 
in string inspired theories, while the third one looks good in grand unification 
models especially keeping in mind that the fine structure coupling constant a 
becomes quite large, about 1/40, at the unification scale. The latter makes the 
mechanism much more efficient. 



4 Generation of vorticity perturbations by cosmological inhomo- 
geneities 

There are several mechanisms which might create inhomogeneities in the primeval 
plasma whose relaxation could lead to generation of primordial magnetic fields: 

1. First order phase transitions creating bubbles of one phase inside an- 
other Ml. Though the magnitude of magnetic fields produced on the 
boundaries between the phases could be very large, the characteristic 
scale is extremely small and it is difficult to stretch it up to galactic size. 

2. Creation of stochastic inhomogeneities in cosmological charge asymme- 
try, either electric EJ ) r e.g. leptonicEJl at large scales which produce 
turbulent electric currents and, in turn, magnetic fields. The first of these 
models could be quite efficient but it demands rather unusual physics. For 
realization of the second model only one rather innocent assumption is 
necessary, namely, an existence of sterile neutrino, v s , very weakly mixed 
with active ones, v e ,p,,T- If v s is lighter than v a and the mass difference 
is sufficiently large then the MSW-resonance transition between v s and 
active neutrinos would take place giving rise to a large, about 0.1, and 
strongly fluctuating lepton asymmetry in the active neutrino sector D3 . 
When the wave length of the domain with a large lepton (or anti-lepton) 
number crossed horizon, the neutrino or antineutrino flux from this do- 
main would induce electric current by scattering on electrons or positrons 
because of different ve~- and i/e + -cross-sections. In this process hydro- 
dynamical flows with large Reynolds numbers would be generated and 
turbulent vortical currents could be produced. The characteristic wave 
length of the generated magnetic field is about 100 pc and chaotic line 
reconnection and large but not unreasonable dynamo amplification are 
necessary for explanation of the observed galactic fields. 

3. Generation of seed magnetic fields slightly before hydrogen recombination 
epoch EJ through relaxation of the usual density perturbations which are 
known from observations to exist. This model does not demand any new 
physics and predicts quite promising amplitude of seed fields at galactic 
scales. Since the work EI) 

is new, not yet published, we will discuss it in 
the next section in some detail. 



5 Generation of large scale magnetic fields at recombination epoch 

We will consider the period when the universe was already quite old and cool 
with temperature about 1-100 eV, i.e. somewhat before hydrogen recombina- 
tion. The usual cosmological density perturbations are known to exist at that 



time, with rather small amplitude, Sp/p ~ 10 4 . The motion of the cosmo- 
logical plasma, or better to say fluid, under pressure forces is governed by the 
hydrodynamical equation (see e.g. the bookE3>) : 



p (d t Vi + v k d k Vi) = -dip + d k 



V d k Vi + div k - -8 ik djVj + di (( djVj) 



(7) 



where v is the velocity of the fluid clement, p and p are respectively the energy 
and pressure densities of the fluid, and r\ and C, are the first and second viscosity 
coefficients. In the case of constant viscosity coefficients this equation is reduced 
to the well known Navier-Stokes equation. The coefficient r\ is related to the 
mean free path of particles in fluid as 



v/p = v = lf 



(8) 



In what follows we disregard the second viscosity £. 

The behavior of the solution to eq. 10 crucially depends upon the value 
of the Reynolds number 



R\ 



v 



(9) 



where A is the wavelength of the velocity perturbations. If R » 1, then the 
fluid motion would become turbulent and non-zero vorticity would be created 
by spontaneously generated turbulent eddies. In the opposite case of low R the 
motion is smooth and and no vorticity could be generated in the first order in 
5 p. One would expect that in the case of scalar perturbations vorticity remains 
zero in any order in 5p. However, this is not the case, as we will argue in what 
follows. 

Let us first estimate the Reynolds number of the fluid motion created 
by the pressure gradient in cq. 0). To this end we assume that the liquid is 
quasi incompressible and homogeneous, so that the second term in the r.h.s. of 
this equation can be neglected. This is approximately correct and the obtained 
magnitude of the fluid velocity is sufficiently accurate. In this approximation 
eq. 0) reduces to a much simpler one: 



<9(V + (v V) v — t/Av 



Vp 
P 



(10) 



A comment worth making at this stage. The complete system of equa- 
tions includes also continuity equation which connects the time variation of 
energy density with the hydrodynamical flux (see below eq. (|19|l ) and the Pois- 
son equation for gravitational potential induced by density inhomogeneities. 
We will however neglect the gravitational force and the back reaction of the 



fluid motion on the density perturbation. This approximation would give a 
reasonable estimate of the fluid velocity for the time intervals when acoustic 
oscillations are not yet developed, i.e. for t < X/v s , where A is the wave length 
of the perturbation and v s is the speed of sound (in the case under considera- 
tion v 2 = 1/3). In fact the wave length should be larger than the photon mean 
free path, to avoid diffusion damping, and the characteristic time interval, as 
we see below, should be somewhat larger than A. So we may hope that our es- 
timates of the velocity are reasonable enough. Neglecting gravitational forces, 
especially those induced by dark matter would result in a smaller magnitude 
of the fluid velocity, so the real effect should be somewhat larger. 

For small velocities (or sufficiently small wavelengths) we may neglect the 
second term in the l.h.s. with respect to the third one. In this approximation 
the equation becomes linear and can be easily solved for the Fourier transformed 
quantities. Assuming that the parameters are time-independent (though it is 
not necessary) we obtain: 

v k = - J^ 5 k [1 - exp(-^ 2 t)] (11) 

where 8k = (5p/p)k is the Fourier transform of relative density perturbations, 
Sp/p; its natural value is ~ 10~ 4 , though it might be much larger at small 
scales. The coefficient 1/3 comes from equation of state of relativistic gas, 
V = P/3- 

Therefore, for the Reynolds number we obtain: 

R * = T7T^ - cxp (- vk2t )] ■ ( 12 ) 

3 (fez/) 

If Sk is weakly dependent on k, then Rk is a monotonically rising function of 
the wavelength \ — 2ir/k. For t ^ \ 2 /v it takes the value 

flfc = ^-5 fc <l, (13) 



so 

gener 



the hydrodynamical flow remains laminar and vorticity is not spontaneously 
aerated. Dynamical generation of vorticity is governed by the equation: 

d t n -i/An = -vx ( — ) (w) 



p 

where fi = V x v and we assume that velocity is small so that the term 
quadratic in v was neglected. If the r.h.s. is non-vanishing, then CI would be 
non-zero too. However usually pressure density is proportional to the energy 
density, p = wp, with a constant coefficient w and hence V x (Vp/p) = O.We 



can see that this is not so because cosmic plasma consist of different compo- 
nents whose motion is somewhat different. Let us assume that plasma is in 
local thermal equilibrium with common temperature r(x). This assumption is 
justified by a large interaction rate between radiation and charged particles. If 
T would be the only parameter which determines the state of the medium, then 
vorticity would not be generated because we would have in our disposal only 
VT and it is impossible to construct non-vanishing V x v from the gradient 
of only one scalar function. However, distributions of charged particles depend 
upon one more function, their chemical potential: 



/ = exp 



T(x 



(15) 



where the dimensionless chemical potential £ can be readily expressed through 
particle number density n e ~ tlb = /3(x)n 7 with (3(x) = 6 • 10~ 10 + Sf3(x): 

£(x) = In fi(x) + const (16) 

Hence we will find that the source term in the vorticity equation (|14|) is equal 
to: 

S k = -e^k 9j\ — = eij k - — - (17) 

VP/ 3 /°tot P Ptot 

An essential feature here is that the spatial distribution of charged particles 
does not repeat the distribution of photons and hence the vectors Vp 7 and 
V/3 are not collinear. This could occur if, for the wave length corresponding 
to subgalactic scales, there exist baryon isocurvature fluctuations and thus 
p(x) and (3(x) have different profiles. As we have mentioned above, different 
mean free paths of photons and charged particles would maintain such non- 
collinearity of the order of unity at the scales A ~ Z 7 . Moreover, even in the case 
of adiabatic perturbations a shift in the distribution of photons and charged 
particles could also be created because of acoustic oscillations that proceeded 
with different phases of radiation and matter densities. At the scales A < Z 7 
perturbations in the the plasma temperature would be erased by the diffusion 
damping , while for A 3> Z 7 the diffusion processes are not efficient and one 
would expect self-similar perturbation leading to collinearity of Vp 7 and V/?. 
On the other hand, when A entered under horizon acoustic oscillations begun 
which destroyed the self-similarity. Thus the expected wavelengths of vorticity 
perturbations should be between l 1 < A < if -1 . 

Surprisingly vorticity can be also generated (and in the case under consid- 
eration even a larger one) if perturbations in plasma are determined by a single 
scalar function, for example, by T(t, x) because it might be proportional to the 
product diT(t,x) djT(t' , x). These two gradients generally are not collinear if 



taken at different time moments t and t' . To see that, let us start from the 
Boltzmann equation for the distribution function f(t,x,E,j>) of photons: 



-+V-V-Fp— + P— J f(t, x, E, p) = I coll [/„, / 6 , ...] , (18) 

where V = p/£7 is the particle velocity (not to be confused with the velocity 
v of macroscopic motion of the medium, for photons V = 1, while v <C 1), 
E and p are respectively the particle energy and spatial momentum, H is the 
universe expansion rate, F is an external force acting on particles in question 
(the latter is assumed to be absent), and / [f a , ft,, ...] is the collision integral 
depending on the distributions f a of all participating particles. 

At temperatures in eV-range only the Thomson scattering of photons on 
electrons is essential, so the collision integral is dominated by the elastic term. 
Integrating both parts of eq. (|18|l over <i 3 p/(27r) 3 we arrive to the continuity 
equation: 

?i(x) + VJ = (19) 
where J is the photon flux given by 

and v is the average macroscopic velocity of the photon plasma. Using the 
standard arguments one can derive from eq. (|19|l the diffusion equation: 

n = DAn (21) 

where D k Z 7 /3 is the diffusion coefficient. We will use this equation below to 
determine time evolution of the photon temperature T. 

If the elastic reaction rate r e i = aThn e X e — 1/Z 7 is sufficiently large, local 
thermal equilibrium would be established and the photon distribution would 
be approximately given by 

/«/o = exp(-£/T + £) (22) 

where the temperature and effective chemical potential could be functions of 
time and space coordinates: T = T(t,x) and £ = £(£,x), and the photon 
mean free path is given by L f = 30 pc /X e (T)T^ v , where T e y is the plasma 
temperature in eV and X e is a fraction of the free electrons: X e (z) is practically 
1 for z > 1500, and sharply decreases for smaller z's, reaching values ~ 10~ 5 
at z < 1000. 
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Evidently /o annihilates the collision integral. We can find correction to 
this distribution, / = /0 + /1, substituting this expression into kinetic equation 
l|18|) and approximating the collision integral in the usual way as — r e i/i: 

(K + T cl ) fx = -Kf (23) 

where K is the differential operator, K = dt + (VV). The solution of this 
equation is straightforward: 

r r T i 

2 T e i (t - t 2 ,x - Vt 2 ) 

^(i-n.x-Vn) (24) 

Using this result we can calculate the average macroscopic velocity of the 
plasma. The calculations are especially simple if elastic scattering rate is high 
and the integrals are dominated by small values of t\. In this case we obtain: 

V ' {t,) ~ f<Ppf (t,x,E) (25) 
and the vorticity, fli = CijkdjVk is: 

To estimate time derivatives of the temperature we will use the diffusion 
equation (|21|l . from which we find d{F = DAT and finally obtain for vorticity 
with the wave vector k = 2ir/\: 

2 



(f ) t»fc*«3.10"(^ ) ^ (27) 



Since the photon diffusion erases temperature fluctuations at the scales A < l~ f 
the vorticity reaches maximum value near A ~ Z 7 . This magnitude of vorticity 
is considerably larger than found previously with the source term (|17f) and we 
will rely on it in the estimates of magnetic field presented below. 
Since the conductivity of cosmic plasma is very high, 

K = (3/2a) K/n 7 ) (m 2 e T), (28) 

the generation of magnetic field by the source currents, created by the cos- 
mological inhomogeneities, is governed by the well know equation of magnetic 
hydrodynamics : 



d t B = Vx(vxB) + -VxJ 

K 



(29) 



The electric current J induced by the relaxation of the density inhomogeneitics 
would contain two components: electronic and protonic. However the first one 
is surely dominant because it is much easier to drift electrons than heavier 
protons. This is why a non-zero current can be induced in electrically neutral 
medium. Of course the motion of electrons would not produce any excess of 
electric charge because the current could be realized by the flow of the dominant 
homogeneous part of electron distribution. 

The solution of eq. I|29|l can be roughly written as 

J""" < 30 > 

An estimate of magnetic field without pregalactic dynamo enhancement can 
be easily done if the helical source current is known, V x J = en e fl. With ft 
given by eq. (123 and B by eq. H31 we obtain: 

*^0^a0»(W(i)(t) 3 (£) 2 ,10-^ (31) 

where we took the wavelength equal to the photon mean free path, A = Z 7 . 

If we take into account that linear compression of pregalactic medium in 
the process of galaxy formation is approximately r ~ 10 2 , the seed field in a 
galaxy after its formation would be r 2 Bo, i.e. 4 orders of magnitude larger than 
that given by eq. I|31|l and, for T = 1 eV, a relatively mild galactic dynamo, 
about 10 4 , is necessary to obtain the observed galactic magnetic field of a few 
micro-Gauss at the scale Ib ~ (100/r) kpc = 1 kpc. The seed magnetic fields 
formed earlier (at higher T) would have larger magnitude (<~ T 3 ) but their 
characteristic scale would be smaller by factor 1/T 2 . Chaotic line reconnection 
could create magnetic field at larger, galactic scale l ga u but the magnitude of 
this field would be suppressed by Brownian motion law - it would drop by the 
factor (Is/lgai) 3 ^ 2 - It is interesting that according to these results all scales 
give comparable contributions at l ga i. This effect may lead to an enhancement 
of the field but it is difficult to evaluate the latter. Let us also note that 
magnetic fields generated by the discussed mechanism at the cluster scale, 10 
Mpc, should be not larger than 10~ 8 p,G if no additional amplification took 
place. 

Larger density perturbations could be helpful for generation of larger 
magnetic field for which dynamo might be unnecessary. Though much bigger 
ST is not formally excluded at the scale about 100 kpc, but to have them at 
the level (ST/T) 2 ~ 10 -4 seems to be too much. A natural idea is to turn to 
a later stage, to onset of structure formation when Sp/p becomes larger than 
10~ 2 . With such density perturbations strong enough magnetic fields might 
be generated without dynamo amplification. However after recombination the 



number density of charge carriers drops roughly by 5 orders of magnitude. 
Correspondingly Z 7 rises by the same amount and the strength of the seed 
field would be 5 orders of magnitude smaller if density perturbations and the 
temperature of formation remained the same. However both became very much 
different. Density perturbations rose as scale-factor, (Sp/p) 2 ~ (T eq /T) 2 , where 
T eq ~ 1 eV is the temperature when radiation domination changed into matter 
domination and density perturbations started to rise. Since, B/T 2 ~ T 3 , 
according to eq. (|31fl . the net effect of going to smaller T is a decrease of B/T 2 
which would be difficult to cure even by later reionization. Still, as argued in 
ref. , magnetic field generation, driven by anisotropic and inhomogeneous 
radiation pressure (and in this sense similar to our mechanism) at the epoch of 
reionization, could end up with the field of about 8 • 10~ 6 pG. This result is 8 
orders of magnitude larger than that found in the earlier papers and quite 
close to ours (JSTJ; though these two mechanisms operated during very different 
cosmological epochs and were based on different physical phenomena. 

Generation of magnetic field at recombination was also considered in 
re f 119J) wnere a much weaker result was found. This difference can be possibly 
attributed to the following effects. We considered above an earlier period when 
the photon mean free path was much smaller than the horizon. It gives a factor 
about 10 3 in fluid velocity, ea. l|ll|) . Moreover, since in our case the electrons 
are tightly bound to photons the electron-photon fluid moves as a whole (while 
protons and ions are at rest) and the electric current induced by macroscopic 
motion/oscillations of plasma is noticeably larger. 

6 Conclusion 

Despite many suggested models, the origin of galactic magnetic fields and, 
especially, intergalactic, if existence of the latter is confirmed, remains myste- 
rious. One class of models is based on known physics and do not invoke any 
ad hoc assumptions for the explanation of the phenomenon. The explanation 
based on stellar ejecta possibly encounters serious difficulties because of energy 
constraints. The mechanism of field generation just before hydrogen recombi- 
nation looks reasonably good but probably it cannot explain both galactic and 
intergalactic fields within the frameworks of the standard cosmology with flat 
spectrum of density perturbations. 

Another class of models is based on physical phenomena in the early 
universe and its different members are spread between inflationary stage to 
relatively late MeV-epoch. These models manipulate with unknown physics 
(except possibly the MeV one) and because of that are much less restricted in 
their possibilities. It is a difficult task to understand what mechanism is indeed 
responsible for creation of the observed magnetic fields. A critical test would 
be a possibility of simultaneous explanation of galactic and intergalactic fields. 



To this end a confirmation of a possible existence of the latter is of primary 
importance. 
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